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Abstract 

We study theoretically the entropic elasticity of a semi-flexible polymer, such as DNA, confined 
to two dimensions. Using the worm-like-chain model we obtain an exact analytical expression for 
the partition function of the polymer pulled at one end with a constant force. The force-extension 
relation for the polymer is computed in the long chain limit in terms of Mathieu characteristic 
functions. We also present applications to the interaction between a semi-flexible polymer and a 
nematic field, and derive the nematic order parameter and average extension of the polymer in a 
strong field. 
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I. INTRODUCTION 



Mechanical properties of biomolecules, such as their response to an applied force, are 
important for understanding a variety of biological processes ranging from cell motility to 
gene regulation. Detailed mechanical studies of biopolymers using laser tweezers, magnetic 
tweezers, atomic force microscopy, and other single molecule techniques have provided the 
necessary experimental input for formulating precise mathematical models of their elasticity. 
In the case of double-stranded DNA (dsDNA) the favored theoretical model is the worm-like 

n 

chain (WLC), proposed by Kratky and Porod in 1949 pj. Owing to the central role that 
DNA plays in biology, the last decade has seen a significant body of theoretical work on 
the WLC model, as well as its modifications and extensions. 

Considerable experimental evidence has now accumulated that shows that the worm-like- 
chain model reproduces the mechanical behavior of dsDNA in the entropic regime quite 



well 



Q,y, 



The effect of confinement on the statistical properties of polymers is of growing interest 
to scientists and engineers [12I . In the context of the cell, the important observation 
is that dsDNA is always found in a state of confinement. Namely, it occupies a volume 
which is considerably smaller than the volume it assumes free in solution. Here we examine 
dsDNA confined to two-dimensions. While not immediately of relevance to cell biology, this 
type of confinement has been studied recently as an interesting polymer physics problem. 
Experiments have probed dynamics and thermodynamics of dsDNA confined to a mica sur- 
face [blj] and the surface of a lipid bilayer 



15] . More generally, two-dimensional confinement 



is emerging as an important experimental method, as it allows for the use of fluorescence 
microscopy and atomic force microscopy to obtain images of individual dsDNA molecules. 
For example, recent experiments have made use of dsDNA absorbed on treated mica surfaces 
to study the effect of DNA-binding proteins on dsDNA conformations [3] . 

In this paper we study the effect of two-dimensional confinement on the entropic elasticity 
of a worm-like-chain polymer. In particular, we compute the partition function for the WLC 
model in the presence of an applied force when the chain is restricted to two dimensions. 
This allows us to derive an exact expression for the average end-to-end distance of the 
polymer as a function of applied force, which is the central result of this paper. 

The paper is organized as follows. In sections [Til and II III we introduce the WLC model 



and derive an exact closed-form expression for the tangent partition function and the force- 
extension relation of a WLC polymer pulled at one end by a constant force, in the limit 
of large polymer length. In section IIVI we discuss the asymptotic form of our result in the 
limit of strong and weak force, and compare it to the approximate force-extension relation 
in three-dimensions. We also give a simple algebraic approximation to the force-extension 
relation in two-dimensions. In section |V] we show that a minor transformation converts the 
partition function of section [Til into the partition function of a semi-flexible polymer in a 
nematic field. We use this to calculate the nematic order parameter for the polymer, as well 
as its relative extension in a strong nematic field. 

The results reported here should be of practical use in stretching experiments where semi- 
flexible polymers such as DNA are confined to a surface. On the purely theoretical side, 
they establish an interesting connection between the elastic properties of two-dimensional 
semi-flexible polymers and Mathieu functions. 

II. THE WORM-LIKE CHAIN MODEL 

The worm-like chain model describes the polymer as an elastic filament characterized by 
a rigidity parameter k, which has dimensions of length. For a three-dimensional polymer k is 
the persistence length, i.e., the decay length associated with the tangent-tangent correlation 
function, which for DNA under physiological conditions is about 50nm. The energy of 
bending a small arc of length s of the filament into a circular segment of radius R is given 
by kBTns/2R 2 , where fcgT is the thermal energy (at room temperature ksT = 4pNnm). 
If we describe polymer conformations by a smooth continuous curve, the Hamiltonian of 
the WLC can be expressed as the integral of this segmental bending energy over the entire 
polymer curve. Therefore, when a polymer is tethered at one end and pulled with a constant 
force F in the cc-direction at the other, the Hamiltonian H can be written as: 

where t is the unit tangent vector to the polymer curve, and L is the contour length of the 
polymer. In writing eqn. (Op) we have made use of the standard relation between the radius 
of curvature of a curve and the rate of change of its tangent vector, \dt/ds\ = 1/R. 

Since this polymer is represented by a curve, its partition function is the sum of Boltzmann 



weights for all possible curves, subject to the constraint of fixed polymer length L. In terms 
of the tangent vector t(s) the partition function can be written as a path integral, 



Z{tf,U-L) = / V[t(s)}exp 




-ft -A 



(2) 



where / = F ■ x/ksT is the reduced force with units of inverse length. The dimensionless 
combination nf delineates regimes of high (nf ^> 1) and low force (/c/ < 1), which we take 
up in section HVl Experiments on stretching polymers that involve attaching the two ends of 
the polymer molecule to beads, constrain the first and final tangent vectors, usually making 
them lie along the direction of the force. This is indicated by specifying the initial and final 
tangent vectors, and t,, explicitly in eqn.(T2]). In case the ends are free, we will need to 
integrate over these tangent vectors to obtain the appropriate partition function. This is 
the case, for example, when the polymer is dissolved in a nematic solvent, as we shall see in 
section IVl 

From the partition function, eqn.(j2J), we calculate the average end-to-end extension of 
the polymer in the direction of the force, (X). Namely, (X) is the conformational average 
of the end-to-end vector R = J Q L dst projected in the direction of the applied force, i.e., 



(X) = ( / dst-x). 
'o 



From eqn.(j2]) we conclude that 



(X) 



d\nZ 



(3) 



(4) 



In two dimensions i = (cos 6, sin 6), where 6(s) is the polar angle, and the partition 
function Z in eqn. (J2J) can be rewritten as, 



Z(6 fl 6 l ;L)= / V[6(s)}exp 



ds I — ( — | — fcosO 
o I 2 \dsj 



(5) 



where 6f and 9i are now the final and initial tangent angles. The calculation of Z(9f,9f,L) 
is made simple by the standard connection between the path integral and the Schrodinger 
differential equation [17] 







1 d 2 



ds 2k 89 2 
subject to the initial condition 



fcos9) Z(9,9i;s) = (s > 0), 



(6) 



Z{9,9 i ;0) = 8(9-9,). 



(7) 



Since eqn.([n]) is separable in s and 9, we can express its solutions in terms of the eigenfunc- 
tions Q{6) and eigenvalues, E, of the corresponding eigenvalue equation, 

1 d 2 



2k dO 2 



+ fcos6)Q = -EQ. (8) 



Note that the partition function, Z(9,9i;s), is formally equivalent to the path integral 
for a quantum rotor with moment of inertia K in an external field /, evolving in imaginary 



18|. The 



time s [2[, or, equivalently, to a quantum pendulum in a gravitational field 
correspondence is actually very intuitive, since thermal fluctuations of the tangent vector 
along the backbone of the polymer can be regarded as the time evolution of the quantum 
rotor in the external field. Each polymer conformation hence represents one particular time 
evolution of the quantum problem, in other words 'the path is the polymer' The 
procedure for computing Z(9,9f,s) we have outlined above is thus formally equivalent to 
obtaining the spectral representation of a quantum propagator in terms of the eigenstates 
of the associated Hamiltonian. 



III. FORCE-EXTENSION IN TWO DIMENSIONS 



The eigenvalue equation, eqn.flE]), can be transformed into the canonical Mathieu differ- 
ential equation by the change of variables @ = 7r/2 — 6/2. As 9 goes from tt to — tt about 
the direction of the force, varies between and tt. This transformation yields, 

d 2 Q 



+ (p- 2gcos26)Q = 0, 



(9) 



with 



8Ek and q = 4/k. The solutions of this equation are the Mathieu functions 



2CJ, |2l|]. The Mathieu differential equation always has a periodic solution and an aperiodic 
solution [2^]. Since rotation of tangent vectors by 2tt must leave Q unchanged, we desire a 
solution that is periodic in 9 with a period of 2tt, or, equivalently, periodic in with period 



TT. 



Solutions of eqn.([H]) with a period of either tt or 2tt are called "basically periodic" in 
the Mathieu function literature. A Mathieu function is basically periodic only for specific 



values of the parameters p and q 20|, 



21 



22]. The value of p that makes the function 



basically periodic is called the "characteristic value" of the Mathieu function, and is in 



fact a continuous function of q |20|, 



21 



22( | . For a given q there is a countable set of such 



characteristic-value functions, usually denoted as a n (q) and b n (q), n = 0,1,2,.... The 
functions a n (q) are associated with the even Mathieu functions, ce n (a n ,q, O), while the 
functions b n (q) are associated with the odd Mathieu functions, se n (b n , q, 6); ce and se stand 
for cosine elliptic and sine elliptic, respectively [20 
Mathieu characteristic functions satisfy, 221 ] . 



21 



22]. Furthermore, for g ^ 0, the 



ao(q) < h(q) < ai(q) < b 2 (q) < a 2 (q) < 



(10) 



We conclude that the periodicity of Q(Q) imposes the quantization condition on the param- 
eter p = 8Ek, which must be from the countable set of Mathieu characteristic functions. 

The periodicity of the Mathieu functions for different values of n can be deduced from 
their expansions in sine and cosine series. There are four types of basically periodic ce and 
se functions, and their expansions are 

oo 



ce 2n (a 2n ,q, 6) 
ce 2n+ i(a 2n+ i, q, 0) 
se 2n+ i(b 2n+ i,q,Q) 
se 2n+2 (b 2n+2} q,Q) 



1 2;- 



cos 2r6, 



r=0 
oo 



E4f + t 1) cos(2r + l)e, 



r=0 



r=0 
oo 



Ei^ 2) sin(2r + 2)e. 



(11) 
(12) 

(13) 

(14) 



r=0 



Note that the coefficients A and B are functions of q. Perusal of the above equations shows 
that the class of functions that have period tc in 0, are the ce 2n and the se 2n+2 functions. 



The normalization convention for the Mathieu functions is 



2Q|, 



ce 2 (a„, q, Q)d 9 



se 2 (b n , q, Q)d 6 



7T 



(15) 



The first few basically periodic Mathieu functions are plotted in fig.flT]), while the corre- 
sponding characteristic functions, are plotted in fig.©. 

With the eigenfunctions [ce 2n and se 2n+2 ) and the corresponding eigenvalues (a 2n and 
b 2n +2) of eqn.(jHJ) in hand, the spectral representation of the partition function of the worm- 
like chain is 



n=0 



Z(Qi,Q f ,L) = ^2<-ce 2n (a 2n ,AfK,Qi)ce 2n (a 2ni AfK,Q f )exp 



-— a 2 „(4//«) 



+ -se 2n+2 (b 2n+2l 4/ac, ®i)se 2n+2 (b 2n+2l 4/k, 9/) exp 

7T 



-—b 2n+2 (AfK) 



16) 




Angle [radians] 



FIG. 1: (Color online) Eigenfunctions ceo, se 2 and ce 2 plotted for 4/k = 1 and for between —tt 
and 7r. 

where Bj and 0/ specify the tangent vectors at the two ends of the chain. To obtain the 
partition function with no constraints on the initial and final tangent vectors, we integrate 
over all possible values of the initial and final angles. Using the expansions, eqns. (fnj ).(TU 
and (113]) . we get 



£(L) = J> (A? ) exp 



n=0 



-^-Ci2n(4/ft) 



(17) 



From the inequalities in eqn.t flOl) it follows that, for L large compared to n, the above 
expansions for Z are dominated by the term containing ao(4fn) in the exponential. Even 
for relatively short polymers, this term is a good approximation to the full partition function. 
Namely, as can be seen from fig.(T5]), the difference between 6 2 (4/k) or a 2 (4/«:) and a (4fK) 
is 4 at f — and increases in absolute value with increasing /. Therefore, even for polymer 
length L w 8k, the second exponential in the expansion in eqn. (1161) and eqn. (U7l) is at least 



FIG. 2: Characteristic functions ao(q), &2(<?), 02(9), b^q) and 04(5) plotted against q. 



e -4 « 0.02 times smaller than the first; the subsequent terms are exponentially smaller. 

It remains to consider the coefficients. In the case of eqn. ( 1T6]) the coefficients are products 
of Mathieu functions and when df ~ 9j ~ 0, which is a reasonable assumption for most 
experiments, the coefficient of the first term is larger than any of the subsequent ones. In 
the case of eqn. f|T7|) . the coefficients (A^) 2 are complicated functions of the dimensionless 
force q. However, all of them satisfy (A^) 2 < 1 and they are all of the same order of 
magnitude when q > 1. When q = 0, A^ = 1/V2, and all other coefficients are zero 20]. 
We conclude that the approximation, 

Z « cexp[-^o (4/K)], (18) 

OK 



where c is the appropriate coefficient from either eqn. ([16]) or eqn. ffTTj) . depending on the 
boundary conditions, is justified for polymer lengths L > 8k. 

We can now calculate the relative extension from eqn.(|3j). The coefficient of the expo- 
nential yields a correction which is of order k/L and hence can be safely ignored in the long 



chain limit. The resulting force extension relation is therefore given by, 

(X) 1 da {4:fK) 1 da {q) 



(19) 



df 2 dq ' 

where q = 4/k, as before. 

The Mathieu Characteristic functions are interesting in themselves. So far no exact ex- 
pression for them is known in terms of other functions, though good polynomial approxima- 



tions are available for small and large 



q- 



involving infinite continued fractions 20 



he y a re usually expressed as roots of an equation 



21 



22 



. They can also be shown to be the roots of 



an infinite tri-diagonal determinant equation [2l|. Both these representations allow numer- 
ical calculation as well as numerical differentiation by mathematical programs like Mathe- 
matica 23|. Mathematica executes the desired function as MathieuCharacteristicA [r ,q] , 



which we use to generate the plot of the relative extension as a function of the applied force, 
shown in fig.Q. 

It is worthwhile to note that the higher order moments of the probability distribution 
of the end-to-end distance can be obtained by successive differentiation of the partition 
function, 

1 r) n 7 

(* n ) = (20) 



ZdF n ' 

Having computed the exact expression for the force-extension relation for a two- 
dimensional worm-like chain in terms of Mathieu functions, in the next section we consider 
the small and large force limits of eqn. (fl~9]) . which are readily obtained using more elementary 
methods. 



IV. LIMITS OF SMALL AND LARGE FORCE 

In order to gain a better understanding of the derived force-extension relation, eqn. ffTQl . 
we analyze it in the small force (q = 4/k <C 1) and the large force (q ^> 1) limits. 



For small forces, we use the power series expansion 22] of ao(q) 



«o(g) = -^ + ^ + 0(g 6 ). (21) 



Therefore, to first order in q, the extension becomes 

(X) 1 



L 2 



q = 2fK. (22) 




1 



2-d approx. 

2-d exact 

3-d approx. 
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FIG. 3: Force extension curves for semi-flexible polymers. The full line is the relative extension 
(X)/L plotted against the dimensionless force, q = AFK/ksT using the exact formula, eqn. (fT9]) . 
The dashed line is the interpolation formula, eqn. (|33p . For comparison, the force-extension curve 
in three-dimensions is plotted as a dot-dashed line, using the approximate interpolation formula 
described in Ref. 



It is reassuring to find that a perturbation analysis of eqn.(jH|), in the limit of low force, 
yields the same result. If the potential / cos 9 is treated as a perturbation, the unperturbed 
equation describes a free particle on a ring, 

1 d 2 Q n 



2k d9 2 



(23) 



The solutions of this equation are sine and cosine functions, except for the = solution 



(24) 



which is Qo = l/y2ir. The eigenvalues are given by, 

,2 



n 
2^? 



n = 0,1,2, 



The first order correction to the zero-eigenvalue is zero. Using orthogonality properties of 
sine and cosine functions, we obtain the second order correction 

E® = -nf. (25) 

Using (X)/L = —dEo/df, which follows from eqn.(|4]), and the fact that in the long-chain 
limit Z ~ exp(— E L), we get eqn. f[2~2l . 

Let us now consider the limit of large forces. We use the asymptotic expansions given in 
Ref. 



221, 



a (q) ~ -2q + 2q^ - 1 - -L g 4 + O^ 1 ). (26) 



This gives us, after ignoring terms of order q 1//2 and smaller, 



(27) 



This formula can be derived independently by noting that at high forces the polymer is 
highly extended in the direction of the force, with only small transverse fluctuations, hence 
the angle 9 that t makes with x is very small [2j. We can therefore use the approximation 

cos0«l-~0 2 . (28) 

This transforms the differential equation to be solved, eqn.flS]), into a form that is equivalent 
to the Schrodinger equation for a one-dimensional simple harmonic oscillator with h = 1 
and angular frequency on = ^/ f / k, 

1 d 2 Q n 1 



2k Ott* " 2 f6 Qn = ~ {K + f)Qn - (29) 
The eigenfunctions Q n are gaussian functions multiplied by Hermite polynomials, and the 
eigenvalues are given by the well known expression, 

E n + f=(n + ^)y/ffa. (30) 
Therefore the partition function can be written as, 

oo 

Z(0 f , 9i\ L)=J2 Qn(9 i )Qn(0 f ) exp [-E n L] . (31) 

n=0 

Taking only the first term, in the limit of large-L, and integrating over all values of the 
initial and final angles, we get, 



Z ~ exp[— EqL] = exp 



(32) 



Using this result in eqn.pi immediately gives eqn. f[2"T|) . as expected. 

In fig. ([3]) we have plotted a comparison of the force-extension curves from eqn. ffl9l) and 
the interpolating formula of Ref . Q ■ The 2d curve is found to lie strictly above the 3d curve 
at non-zero force. This is expected by entropic arguments since a three dimensional polymer 
has larger entropy, hence a larger force is needed to extend it by the same amount. Note 
that in two dimensions DNA is highly stretched even at relatively low forces. 

It is also possible to write down an approximate formula for the force-extension of a two- 
dimensional polymer, which tends to the exact force-extension relation in the two limits, 
and is approximately true to within less than 10% in the remaining regime (See fig. (j3J)), 

16 /« = - 1 + (33) 



L J 

V. SEMIFLEXIBLE POLYMER IN A NEMATIC FIELD 

Another way to stretch a polymer is to put it in an aligning field, such as that produced 
by a nematic liquid crystal. In a recent experiment, semi-flexible polymers were stretched 
by dissolving them in the nematic phase of rod-like /d-viruses 
experienced by the polymer can be written as H ]. 

V ^ 



241 ] . The nematic potential 



r / ds{{i-nf -1/2} , (34) 



k B T Jo 

where the square arises due to the reflection invariance of the unit nematic director n, 
which we assume again is in the cc-direction. The coupling T, between the nematic field 
and the polymer depends upon the order parameter of the nematic liquid crystal in which 
the polymer is embedded. When the liquid crystal undergoes a phase transition from the 
isotropic to the nematic phase, F increases and the polymer stretches. This mean-field form 



of the interaction ignores the fluctuations in the nematic order parameter 26|, 1271 ] . 

With the help of the trigonometric identity, cos 2 9 — 1/2 = (cos2#)/2, the potential in 
eqn. (1341) becomes similar to that of the constant pulling- force case, but with cos# replaced 
by cos 29. Following the treatment detailed in a previous section, the partition function 
can again be expressed in a spectral representation, using the characteristic functions of the 
Mathieu differential equation, which for this case is, 

d 2 

— £ + TkQ cos 29 = -2EkQ. (35) 
o9 2 



As before the ce 2n and se 2n +2 functions are the appropriate eigenfunctions. In the large- 
ly limit, the smallest eigenvalue again dominates the spectral expansion and the partition 
function is well approximated by, 



L I 




2k 





Z = 27r(4 0) ) 2 exp 

The derivative of the free energy with respect to the coupling T is now 

dlnZ rL 



(36) 



S 



{[ ds{cos 2 6- 1/2}), (37) 
Jo 



dT 

which is a measure of the degree of alignment of the polymer with the external field. It can be 
seen that when the polymer lies perfectly along the x direction, then S — L/2, but when the 
polymer conformation is completely random, 5 = 0, since (cos 2 6) = 1/2. Therefore 2S/L is 
the nematic order parameter of the polymer. The partition function eqn. (1361) enables us to 
calculate 5/L as a function of the coupling strength T, 

L ~ ~4^~' (38) 
where now q = Fk/2, and we have ignored terms of order k/L. Using the above result we 
can also obtain the approximate relative extension of the polymer in the direction of the 
nematic director at high coupling T. Due to the reflection invariance of the nematic director, 
the extension along that direction is defined as, 

(| X |) = (| f dscos9 |). (39) 
Jo 

When the coupling between the polymer and the field is large the polymer conformation 
would consist of roughly straight sections which have small fluctuations in the y direction, 
interspersed by hairpin bends. A simple thermodynamic argument, based on comparison of 
the energy cost of a hairpin with the entropy of a hairpin, leads to the conclusion that the 
typical length of the straight segments grows as the exponential of the square root of the 
coupling strength, 

^VttkT) . (40) 
This is similar to the relation discussed in ref. 27J]. For DNA, k w 50nm, therefore L nem w 



50fim when F = 0.05nm _1 , and L nem m 3m when T is an order of magnitude larger. If 
the length of the polymer is smaller than L nem , hairpin bends are ruled out, and cos^ does 
not change sign along the polymer contour. Without loss of generality, we can take it to 



be positive and remove the absolute value in eqn. (1391) . This allows us to write the following 
approximate formula for the extension, 

(X)n(^ dsy/l- sin 2 6) « (jf d s |l - ^ sin 2 flj). (41) 

By arguments similar to those made earlier for the polymer under a constant force, we know 

that eqn. (1361) is a good approximation for the partition function when L/k » 1, and we 

may use it to calculate the right-hand side of eqn. (141 1) . The relative extension of the polymer 

along x when L < L nem and L/k » 1 can therefore be expressed as, 

(X) _ 3 _ 1 da (q) 
L ~4 8 dq ' 1 j 

where g = Tk/2. In the limit q >> 1, the relative extension has the simple form, 1 — 

1/8^/q. These results can be used to estimate V from experiments that measure the relative 

extension. 

It should be noted that the exponential dependence of L nem on the square root of V 
implies that when > 1, L nem could be quite large, hence it would be easy to satisfy both 
L < L nem and L/k » 1. For example it was estimated above that L nem was as much as 
50fim for DNA when Tk ~ 2.5. At this value of T, eqn. (1421) predicts that the molecule would 
have extended to 87% of its contour length. While to our knowledge, experimental estimates 
of r for DNA embedded in a nematic liquid crystal do not exist at present, the experiment 
of ref. 24| measured values up to T pa 800(;um) _1 for wormlike micelles (k px 0.5yum) in a 
solution of fd- viruses. The regime of validity of eqn. (1421) is therefore expected to be quite 
large. These arguments also imply that a dramatic straightening of the polymer should 
take place as the liquid crystal goes through an isotropic-nematic phase transition. This 
indicates that there is considerable scope for using liquid crystals for stretching DNA in the 
laboratory. 

VI. CONCLUSION 

At present several schemes exist to calculate the partition function of the worm-like 
chain models in both two and three dimensions, to any degree of accuracy . With 

the exception of ref. all use eigenfunction expansions to obtain the partition function 
numerically. Ref. [?| uses a novel method to express the partition function in terms of infinite 
continued fractions. 



This multitude of approaches only underlines the rich physics of the worm-like chain 
model, and also hints at interesting mathematical connections. In the present work, the 
use of Mathieu functions to derive an exact partition function yields a simple closed-form 
expression for the free energy of the worm-like chain in the long chain limit. This allowed us 
to calculate the force-extension relation in two dimensions in terms of Mathieu characteristic 
functions. We also discussed the application of this result to the interaction between a semi- 
flexible polymer and a nematic field, and derived the nematic order parameter and average 
extension of the polymer in a strong field. It should be noted that self-avoidance, which 
is more important in two dimensions than in three, has been ignored in our analysis. We 
address this question, as well as stretching by a force that does not remain constant along 
the contour, such as that produced by an electric field, in a future publication j^ . 
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